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Abstract 

A study of 3-body resonances has been performed in the framework of con- 
figuration space Faddeev equations. The importance of keeping a sufficient 
number of terms in the asymptotic expansion of the resonance wave function 
is pointed out. We investigated three neutrons interacting in selected force 
components taken from realistic nn forces. Three-neutron resonance pole tra- 
jectories connected to artificially enhanced nn forces could be found. The 
final pole positions corresponding to the actual force strengths could not be 
found due to the onset of numerical instabilities. The numerically reliable 
results, however, achieved makes it likely, that three-neutron resonance ener- 
gies will have large imaginary parts and therefore physically not interesting. 
A straithforward application of complex scaling to 3-nucleon systems with 

realistic forces could not be controlled numerically. 
PACS numbers: 21.30.+y, 21.45,+v, 24.10.-i, 25.10,+s 
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I. INTRODUCTION 



Three-nucleon (3N) resonances have not yet been firmly established. The situation up 
to 1987 has been collected in |l|. A recent investigation |2| of the process 3 He(ir~, n + )3n 
found no evidence of a resonance state of three neutrons. Earlier investigations ||, Q 
pointed to resonances in the three neutron (3n) and three proton (3p) systems. The double 
charge exchange process on 3 He was also investigated in [KJ criticizing previous work but 
pointing again to a 3n resonance, now around (20 — i 20) MeV. For the nnp system an 
excited state around 7 MeV with a width of T = 0.6 ±0.3 MeV has been suggested in B. A 
variational study based on complex scaling and simplified nucleon-nucleon (NN) forces was 
carried through in with the prediction of a 3n resonance at E = (14 — % 13) MeV for a 
J w = 3/2 + state. Earlier theoretical studies |§, @ might be useful in various aspects, but 
are not conclusive with respect to the actual position of 3n resonances. A virtual state of 3 H 
with the quantum numbers of the ground state has been investigated in |10| by analytically 
continuing the momentum space Faddeev equation and even estimating three-nucleon force 
(3NF) effects. In the framework of the hyperspherical harmonic method and using complex 
scaling with model potentials a subthreshold 3n resonance for J w = l/2~ has been located 



in 



njj. Recently mathematical foundations have been laid on the analytical continuation 



of the three-body Faddeev equations into unphysical energy sheets |T2J and an application 
thereof to the 4 He trimer appeared in ||13|| . 

It is a theoretical challenge to determine the locations of the lowest lying 3N resonances 
on unphysical sheets based on the most modern nuclear forces. Those forces are very suc- 



cessful in describing 3N scattering observables |T4| which leads to the expectation that the 
corresponding 3N Hamiltonian will also locate 3N resonances with a high degree of reliabil- 
ity. Like for the 3N bound state energy, where some theoretical underbinding results and 
3NF's are used to bridge that gap to the data point, also 3N resonance positions might fi- 
nally be fine-tuned in this manner. These future rigorous results would also be an important 
guidance for experiments and would clarify the hitherto controversial situation. 
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For the 3N system the bound state has been studied extensively in the last decades 
resulting in numerically precise and converged solutions of the 3N Faddeev equa- 



tions fllEf , JET], [0, (T^j. Energy eigenvalues and wave functions have been gener- 
ated with most modern nuclear Hamiltonians in momentum and coordinate space. Also 
numerically exact solutions of 3N Faddeev equations for continuum states with the same 
Hamiltonians are now available in momentum space for a wide range of energies Hl4| . The 



first model solutions also in coordinate space appeared and variational approaches ma- 
tured and provide equally accurate solutions pl| . 



3N resonances, however, have not yet been investigated in a comparable rigorous manner 
and based on realistic NN forces to the best of our knowledge. In this article we would like 
to take first steps into such a direction. 

In the two-nucleon (2N) system the exact form of the wave function in configuration 
space is known outside the range of the potential. In the 3N system only asymptotic ex- 



pansions are known |]2zfl , [ p3} | and methods how to apply them with sufficient accuracy, for 
instance in Faddeev equations, have still to be elaborated in that framework. There are 
first very promising results for 3N scattering above the 3N breakup threshold pi though 



up to now only for simple S-wave model forces. To the best of our knowledge 3N resonance 
wave functions, which asymptotically oscillate and grow exponentially, have not yet been 
determined as solutions of Faddeev equations and realistic forces. First attempts in this 



direction have been done in ||24j| , however with very simple forces. 

In this study we investigate 3-body resonances as solutions of configuration space Faddeev 
equations. The largest obstacle to find the resonance positions with sufficient precision is the 
approximate form of the resonance boundary condition. Ab initio it is not known at which 
coordinate values this form is valid with sufficient accuracy. For resonances this problem 
can be overcome in principle using the complex scaling method (CSM) P5| , |26|| . Thereby 
all relative coordinates are "rotated" into the complex plane leading to bound-state like, 
exponentially decreasing resonance wave function. This "rotation" of the coordinates into 
the complex plane works well for problems in atomic physics with potentials of inverse power 



laws. It will be shown in this article that in nuclear physics a straighforward application of 
CSM in the configuration space Faddeev form leads to serious numerical problems when the 
usual realistic NN interactions of Yukawa type and hard core behavior are used. Whether 
feasible modifications thereof can be found remains to be seen. Right now our experience is 
such, that the resonance energies can be gained with sufficient precision only if one does not 
leave the real r-axis and if one guaranties the exact boundary condition as close as possible 
by using a sufficient number of terms in its asymptotic form. 

In section II we explain in detail the problems encountered when one applies the CSM 
to the nuclear resonances using the 2N system as a play ground. The exact knowledge of 
the boundary conditions for this system allows one to check assumptions which can later be 
applied to 3N resonances. The 3N equations and the method for their solution are presented 
in section III. There we discuss results for 3 bosons interacting with simple, S-wave forces 
and more importantly for 3n's interacting via realistic neutron-neutron (nn) interactions. A 
summary is given in section IV. Information on the technical performance can be found in 
the Appendices A and B. 



II. RESONANCES IN THE 2N SYSTEM 

Searching resonances in a system of two-nucleons is a well defined problem and numeri- 
cally well under control ||27|| . In spite of that we discuss it here and present some results due 
to its pedagogical simplicity, which enables us to present and test approaches which might 
be useful for finding resonances for three or more nucleons. In addition we can introduce 
our notation. 

The 2N resonances are associated with the poles of the S-matrix which are embedded in 
the fourth-quadrant of the complex k-plane (Re(k)> 0, Im(k)< 0). They are solutions of 
the time-independent Schrodinger equation without incoming wave and the outgoing wave 
increasing exponentially at infinity. We introduce the coordinate space partial wave basis 
|r, a 2 >= \r, (ls)jt > where 
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< r'\rlm>= 5{r ' j^ YiJf) (1) 

and where the orbital angular momentum 1 and total spin s are coupled to the total angular 
momentum j. The isospin of the two nucleons is denoted by t. Then the resonance wave 
function , 0(T > ) to a given j and t has the form 

^) = £^M <f \a 2 > . (2) 
i r 

It is written in terms of the reduced wave functions <p a2 (r) which fulfill the radial Schrodinger 
equation 



EIW^S - Ula ; 2 +1) ) - ^V; 2Q ,(r)]^(r) = -^(r). (3) 



The resonance conditions are such that ip has to be regular at the origin and purely 
outgoing. For the reduced amplitudes this amounts to 

Q2 (O) = O (4) 

and 

4> a2 (r) <x rh<tll(kr) (5) 

outside the range of the interaction. Here h\ (kr) is the spherical Hankel function. As 
we shall demonstrate below it is mandatory to either use the exact form of the boundary 
condition (0) or to take into account at least the leading orders to find the resonances in 
the 2N system with sufficient accuracy. In Eq.([|) m is the nucleon mass and k 2 = jriE. 
The complex energy E is embedded in the lower half plane of the second energy sheet. The 
boundary conditions (|) and (||) apply also to bound states with real, negative eigenvalues 
E. The matrix elements of the NN interaction are defined in terms of V a2a ' 2 (r) by 

< ra 2 \V\r'a' 2 >= 5{ " "/V ^r) (6) 

A method of solving Eq. (^) is displayed in [J28| . One chooses suitable grid points in r 
and interpolates Q2 (r) in terms of cubic splines. The boundary conditions (§]) and @ can 



be conveniently incorporated either by the absence of certain expansion coefficients in the 
spline expansion or by suitable relations among them. Then choosing appropriate collocation 
points in r one ends up with a homogeneous algebraic set of equations. They can be handled 



by standard methods [ES 



In systems with more than 2 nucleons the boundary conditions are known only via an 
asymptotic expansion and one has therefore to work with approximate boundary conditions. 
The knowledge of the exact resonance solutions in the 2N system allow now to study the 
question pivotal for 3N systems: how precisely does one need to know the boundary condi- 
tions in order to get the resonance energies with sufficient accuracy? To study this problem 
one can use the series representation for the spherical Hankel functions 

,(P M - f + 1 exp(i{x - 

1 { ' ^ m\(l-m)\2™ x m + l K ) 

and write the resonance wave function outside the interaction region in a form which resem- 
bles the asymptotic boundary condition in the 3N system (see Section III, Eq.(f2~l|)) 

1 a m 

(f)(r) = exp{ikr) £ — ^_ (8) 

m=0 \ Kr ) 

with a m being unknown constants independent from k and r. Here, contrary to the 3N 
system, the number of terms in the series is restricted by the orbital angular momenta 
contributing to the given total angular momentum j of a resonance. 

In Table | we present 2N resonance positions in different NN states where the exact 
boundary condition have been imposed at different matching radii outside the interaction 
region. The NN interaction Reid93 |29| has been used and has been multiplied with a factor 



Xv (shown in column 2 of Table |). Using the exact boundary condition the resonance 
positions do not depend of course on the value of the matching radius. However, when the 
most crudely approximated boundary condition is applied, keeping in Eq. ([8]) only the term 
with m = 0, the approximate resonance positions vary with the matching radius and up to 
20 fm the results are highly insufficient. As seen from Table | it is sufficient to get correct 
energy eigenvalues if one keeps only the two leading terms in the boundary condition (J8[) ( 
for 1=1 this is the exact form). 

6 



For resonances the problem of a possibly approximate nature of the boundary condition 
can be overcome in principle using the complex scaling method (CSM) [p5|1 , ||26|| . This will 
now be illustrated for 2 nucleons because of its simplicity. In the 2N system one replaces 
the relative coordinate r by 

r _, r > = r e ™. (9) 

For real positive angles a the asymptotic form e lkp with Im{k) < can be changed to 
an exponentially decreasing function. Therefore a resonance wave function analytically 
continued to those complex coordinates can be treated like a bound state with its well known 
boundary condition. Then it is probably enough to restrict the boundary condition (H) to 
the first term m = 0. This "rotation" of the coordinate into the complex plane is sufficient 
for problems in atomic physics with Coulomb potentials. In nuclear physics the potentials 
decrease exponentially and the transformation (|9|) turns then T into e -p-( cos ( a )+ lsm ( a )) r _ 
Therefore for a > the range of the potential will increase, which is counterproductive. 
This can be avoided by generalizing ([|) to complex a's. Then 

e-t* r -> exp{-/ie- /m(Q) [cos(Re(a)) + isin(Re(a))}r} (10) 

which for Im(a) < increases the mass ficos(Re(a)) in the exponent to 
and thus choosing Im(a) suitably decreases the range of the interaction again. At the same 
time also the wave function decreases faster due to the factor e~ /m(Q ) in the exponent. Let 
us call that transformation extended CSM (ECSM). 

We firstly illustrate the power of the ECSM in the 2N system using the Reid 93 NN 
potential in the state 3 Po- As we have seen in Table p], an enhancement factor Ay = 5.5 
leads to a resonance state with a complex energy eigenvalue E = (1.42 — i 0.75) MeV. If one 
uses in Eq.(|8|) the first term with m = only the energy eigenvalue is badly presented (see 
Tables | and |IT|) . The result is improved and we can gain the correct result with the required 
precision using CSM with Im(a) = (see Table |TJ). Note that in Table |TI| we stick to the 
asymptotic form (m = 0) in view of the fact that for 3 bodies we also have only approximate 
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forms. If we now allow for negative /m(«)'s then as can be seen in Table |I| we gain the 
correct eigenvalue even choosing a very small matching radius. 

We also see from Table [II] that CSM for increasing Re(a) gets numerically unstable at a 
certain value of Re(a). In this example it happens around 40°. This behavior is intricately 
connected to the present day parametrizations of NN interactions and puts into question the 
applicability of CSM for finding broad resonances in nuclear systems. The parametrizations 
are typically of Yukawa type (one boson exchanges) like in the Reid93 p9f or Reid soft 



core potential pOfl . The analytical continuation to complex r- values leads to an oscillating 



behavior of the potential, which drastically increase with increasing Re(a), with unstable 
numerics as a consequence. We examplify this behavior in Fig. [I] for the Reid93 potential. 
We checked several of the existing NN potentials and found always this problem ( potentials 



according to Gogny J31J, R.de Tourreil and D.W.L. Sprung p2| , H. Eikemeier and H.H. 



Hackenbroich |33 |, and R. de Tourreil |34j). In case of the AV14 |35] and AV18 |36| potentials 



the situation is even worse due to a singularity for complex r- values in their parametrization 
of the hard core part. Also Gauss-parametrizations are not an exception. In this case the 
form exp(—ar 2 ) limits the complex scaling angle Re(a) to Re(a) < 45°. 

One possibility to generalize the CSM in order to avoid the above problems is the so 
called "smooth-exterior" complex scaling (SECS) f26j. The idea is to perform the complex 
rotation only outside the region of interaction. In this way the analytical continuation of the 
potential into the complex coordinate plane is avoided. In practice this can be achieved for 
instance by performing the following transformation (note the printing errors in |26| leading 
to false expressions (2.1.30) and (2.1.32)) 

r — > r = F(r) — > r e lOL for r — > oo (11) 

with 

, , . . 1, \l + exp(2\(r-r ))]\l + exp(-2\(r-r ))]^ , . 

F(r) = r + \exp(ia) — 1 {r + —In-. ^ —, -) ^U. 12 

K ' J1 4A [l + exp(2X{r + r ))}[l + exp{-2X{r + r ))] f K ' 

This transformation leads to the smooth-exterior-scaled Hamiltonian 
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with 



and 



H = fit \a 2 + v ( F ( r )) + ~73TT^~ + ^TT~ 13 

mj A \r)dr' L m j d {r) dr m r z (r) 



dF 

f(r) = — = 1 + [exp{ia) - l]g(r) 

g (r) = 1 + 0.5[tanh(X(r - r )) - tanh(X(r + r ))] (14) 



— = 0.5A[exp(ia) — l][cos/i™ 2 (A(r — r )) — cos/i~ 2 (A(r + r ))]. (15) 
dr 



The case A = corresponds to usual CSM r' = r e ia [26|. Taking r outside the region 



of the NN interaction and choosing a proper value for A a smooth-exterior-scaling path 
results, which avoids the oscillating changes of the NN potential. We present in Table [T 1 1 
the effectiveness of this approach by aplying it to the l D2 resonance state for the Reid93 NN 
potential with Ay = 5.0. It is enough to keep only the first term in the boundary condition 
(H) to reproduce the correct resonance position when a sufficiently large value of the scaling 
angle a is chosen. Now, however, contrary to CSM, no numerical problems arise when large 
values of a are used. In spite of its success in the 2N system we think that the SECS 
cannot be easily applied in 3N systems. We could not overcome problems arising from the 
permutation operator (see Section III). 

III. THREE-BODY RESONANCES 

Now we investigate 3-body resonances as solutions of configuration space Faddeev equa- 
tions. They are often called Kowalski-Noyes equations and have the form [07j| 



(Hq + V - E)ip = -V P^p (16) 

Here Hq is the kinetic energy, V the two-body force, ip a Faddeev component and P the 
sum of two permutation operators such that the total, properly symmetrized wave function 
is \& = (1 + P)ip. Based on standard Jacobi variables it — — r% , — \{%\ — + ^3)) 
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it is convenient to introduce polar coordinates x = pcosO and y = ^psinO. If we express 
the Faddeev component ip(~x*, in terms of partial wave amplitudes \a > which specify 
angular momenta and isospin quantum numbers, then we obtain 

^,t) = J2^^-<£y\a>, (17) 

a X V 

and Eq.([H]) turns into a coupled set of partial-integrodifferential equations 



(A a + -2 E)0 Q (p, 6) - V a0 (pcos6)<f>p(p, Q) = Y, V a p(pcos6) 
n P P 



X 
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£ f + Kp^e^ip^de' (is) 



with 

Aq= ^ | 10,1a 2 *(* + !) A(A + 1) 
a <9p 2 pdp p 2 dO 2 p 2 cos 2 6 p 2 sin 2 

and 0- = |0- f|, 0+ = ||0-f| - ||. 

The coupling results from the permutation operator P and from possible spin dependen- 
cies in V. The evaluation of the permutation operator matrix elements < xya\P\x'y'a' > 
and the resulting expression for the kernel Kp y (9, 6') is presented in Appendix A. 

Similarily as in the 2-body case the resonance conditions are such that ip(x, y) = has 
to be regular at the origin and purely outgoing in all channels. For the reduced amplitudes 
4>(p, 6) this amounts to [ p8| , |38[ 



0(O,0) = 0(p,O) = 0(p,|) = O (20) 



and for large distances 



ikp i 

ftp, 9) - u(x)e^f + -—t(A(9) + —B{9) + ■ ■ ■) (21) 

(kpp &p 

The first term in Eq. ( pi]) occurs only in channels a, where a two-body bound state u(x) 
exists. The second part describes the behavior in the 3-body breakup channel. The corre- 
sponding momenta are 
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q = Jm(E - E 2 ) 



k = V^E (22) 

where m is the nucleon mass and E 2 the two-body binding energy. The boundary conditions 
(p|) apply also to bound states with real eigenvalues E which are located below all two-body 
thresholds E 2 . For resonances E is complex and the eigenvalues E close to the real axis of 
the physical sheet are of interest. They either lie in second sheets which are accessible to 
the right of ^-values (two-body fragmentation cuts) or in the lower half plane of a second 
sheet, which is accessible through the 3N-breakup cut to the right of E = 0. For a recent 
description of the cut-structure see for instance f24}| . 

An efficient method for solving Eq. (|18D is very well displayed in ||28|| . One chooses 
suitable grid points in p and 9 and interpolates 4>(p, 9) in terms of cubic splines. Similarily 
as in the 2-body case the boundary conditions (|20|) and (|2l|) can be conveniently incorporated 
either by the absence of certain expansion coefficients in the spline expansion or by suitable 
relations among the coefficients. Then choosing appropriate collocation points in p and 9 one 
ends up in a homogeneous algebraic set of equations. For realistic NN forces the dimension 
of the matrix is very large and the tensor product method is a key-technique. We use a 



Lanczos type iteration method j|0| and of course take advantage of the sparceness property 
of the matrices. It is advantegeous to work with further reduced amplitudes F defined as 
4>(p,9) = pe lkp E(p,9), if one treats resonances below the cut with Re E > 0. In Appendix 
B we present the details of the above procedure. 

The main problem in solving Eq. (ITS]) is the approximate nature of the asymptotic form 



(^T|). Ab initio it is not known at which p- value this form is valid with sufficient accuracy. For 



a recent discussion see |38j • Also it is a priori not known how many terms should be included 
in the asymptotic expansion (|2lD for a particular 3N resonance state. One suspects of course 
that including more terms one approaches more closely the exact boundary condition. The 
importance of the terms in Eq. (^) going with increasing powers of depends not only on 
the resonance energy and the matching value for p but also on the magnitudes of the a priori 
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unknown functions A(6), B(9), ■ ■ •. For resonances commonly these problems are expected to 
be overcome using the complex scaling method (CSM) [p5|| . However, due to the restrictions 
imposed on this method by the behavior of the NN forces as discussed in Section II we 
found that this method can be applied in the framework of the configuration space Faddeev 
equations only for rather narrow resonances. Thereby all relative coordinates ( x and y in 
the 3-body case) are rotated into the complex plane. This amounts to p — > p' = pe ia where 
a is restricted to rather small values. An example for realistic nn forces is mentioned at 
the end of this Section. However, if the forces are very simple, CSM can work as will be 
illustrated below. 

One is intended to hope that SECSM which proved to be very efficient in the 2-body 
case will help here, too. However, according to our insights up to now this approach leads 
to problems we could not solve when calculating the permutation operator P and when 
performing the transition from cartesian to polar relative coordinates. 

The only possibility left is to look at solutions of Eq. ( |i8|) without leaving the real r-axis 
and taking the boundary condition (21) with a sufficient number of terms into account in 



order to reach as close as possible the exact boundary condition. In this manner one can 
hope to get convergence for the resulting resonance position. 

We checked our configuration space code by calculating the energy of the bound state 
of 3 if nuclei with the Reid 93 NN potential taking different NN force components into 
account. In Table [IV] we present our numbers together with results of momentum space 
Faddeev calculations The very good agreement between both approaches is clearly 

seen. 

In the first step we considered 3 bosons interacting via Gaussian S-wave forces 

V(r) = -120exp(-r 2 ) + l2exp(-r/3) 2 [MeV]. (23) 

They are such that a barrier is built in, which is an obvious and trivial mechanism gen- 
erating a resonance. That interaction does not support a two-body bound state and in 



the asymptotic form ([21]) the first term is absent. For 3 bosons in the J n = + state 
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and interacting by S-wave forces only the kernel K a p = This model has been used 
in |^2| in applying the stochastic variational method. The resonance energy found there is 
E = (8.60 — i 1.84) MeV. In this case CSM works. For instance for a matching radius of 
p = 10 fm and taking the A- and B-terms in (|21| ) into account we achieved the results shown 
in Table |V|. For Re(a) = 40.0° we are close to the result found in but not close enough. 
A smaller matching radius, however, like p = 5 fm does not work. The correct result can 
however be achieved if in addition we turn on Im(a) < 0. This is demonstrated in Table [V]. 
We illustrate in Figs.||] and [| typical resonance amplitudes F working with real coordinates 
or with complex coordinates. In both cases the matching radius has been chosen as 5 fm. In 
the first case this is not sufficient to achieve the correct complex energy eigenvalue, whereas 
in the second case the amplitude is already sufficiently damped at p = 5 fm and the correct 
energy eigenvalue is accessible at such a small matching radius. 

In this simple case one can study the problem how many terms A, B, ■ ■ ■ in the boundary 
condition (^1|) should be taken into account to get the correct resonance position. Of course it 
will depend on the value of the matching radius p. We present in Table [VT] resonance positions 
obtained with increasing number of terms in the boundary condition (^1|) and different 
matching radii. For the details to include arbitrary number of terms in the asymptotic 
condition we refer to Appendix B. It is seen from Table [VT] that for matching radii > 10 fm 
it is practically sufficient to take only the two leading terms going with A and B into account 
in order to get the correct resonance position. We also see that the restriction to the A term 
alone is totally insufficient. Unfortunately this example is far from reality in nuclear physics. 
It does not contain a short range repulsion and the mechanism for creating a resonance in 
few nucleon systems is different in reality, where it is not caused by a barrier in the potential 
itself. 

Now we investigate the case of three interacting neutrons in states with different total 
angular momenta and parity. We allowed the nn force to act in the states 1 5'o, 3 Pq and 
3 P 2 — 3 F 2 . They are of attractive nature while the repulsive ones were neglected. This 
is only a feasibility study and should provide a first orientation. In the first step the nn 
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forces were chosen as a smooth, local NN force of ref. pty , called in the following Gogny 
potential. It is a rather simple approach to NN phase-shifts through forces of central, tensor 
and spin-orbit type supplemented by a L\i term, all of which of Gaussian form. To initiate 
the search process for the position of the resonance energy we artificially enhanced the force 
components in the three states such that all of them support a two-neutron bound state 
at the same energy. An example is shown in Table [VI I| for a J 71 = 3/2~ state of three 
neutrons. For these enhancement factors also the 3n system is bound. Then we reduce the 



enhancement factors (some intermediate cases are shown in Table [VI I|) , which traces out a 



bound state trajectory. In this case we see that the 3n bound state still exists even when 
all 2n subsystem bound states have disappeared. Then we can neglect the first term in 
the boundary condition (|2T|). Decreasing further the enhancement factor (s) (in this specific 
case only for 3 P 2 — 3 F 2 ) the energy eigenvalue goes into the second sheet and traces out a 
resonance trajectory. The search for the resonance position corresponding to the actual nn 
force is finished when the enhancement factors in all nn force components have reached the 
values Aj = 1.0. As we shall see in the case of three neutrons we face numerical instabilities 
before all X's reach the value 1. We can present the resonance trajectory only partially up to 
the point, where we encounter those instabilities. We expect that they are probably caused 
by the difficulty to fulfill the resonance condition, which requires that the exponentially 
decreasing, purely incoming wave is much smaller than the exponentially increasing purely 
outgoing wave. As we see from Table [VI]] the J 71 = 3/2~ 3n bound state disappears when 



the enhancement factors arrive at (\i So , Xs Po , \3 P2 _3p 2 ) = (1.0,1.0,4.35). In Table |VIII| the 
trajectory of the 3/2~ 3n resonance in the second energy sheet is mapped out by decreasing 
the enhancement factor A3p 2 _3p 2 in step of 0.05 starting from A3p 2 _3p 2 = 4.3. The final 
reliable position reached is for A3p 2 _3p 2 = 3.6 with the value = (5.30 — ?4.67) MeV. 

As can be seen this stable position is only obtained when at least 4 terms in the boundary 
condition ( pT[ ) are taken into account. We checked that our particular distribution of p and 
9 points as well as the choice of the matching radius in p does not influence the resulting 



trajectory. The values given in Table |VIII| are based on a matching radius p = 30.0 fm and 
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typically N p = 30 and Ng = 25 points distributed by scaling factors S p = 1.2 and Sg = 0.9 
(see ref. f2"8"|). It is unfortunate that right now we are not able to reach the final position 



with A3p 2 _3^ 2 = 1. However very likely the width will be quite large. Therefore one cannot 
expect to see a visible enhancement around some 3n cm. energy in a reaction producing 
three neutrons in the state 3/2~. 

Similarily to 3/2~ behave the l/2~ and 3/2 + states based on the Gogny potential. 
With proper enhancement factors Ay the three neutrons are bound together with 2-neutron 
bound states for 1 So, 3 Po and 3 P2 — 3 F^. These 3n states remain bound even when the 
2-neutron bound states disappear with appropriately reduced enhancement factors. They 
disappear into the second energy sheet when the enhancement factors reach (1.0,7.5,4.75) 
and (1.0,5.5,4.8) for J n = l/2~ and 3/2 + , respectively. These enhancements factors are of 
course not unique but we worked with those examples. Following their trajectories before 
numerical instabilities set in results in the complex energy eigenvalues shown for the Gogny 
potential in Table |IX|. For the 1/2" and 3/2 + states, similarily as for 3/2", the final, con- 



verged position is obtained only when at least 4 terms in the boundary condition are taken 
into account. In Fig.^ we show the trajectories for the three states l/2~, 3/2~, and 3/2 + 
for the Gogny potential. We note that the trajectories are not unique and depend how the 
enhancement factors are changed. Our startegy was to bring all of them together to one as 
close as possible. 

In this article we did not investigate the state l/2 + , since we found, like in || that the 
3n bound state disappears earlier than the 2n bound states. In other words, the 3n bound 
state trajectory will enter into the 2n-n two-body fragmentation cut. We left that study 
with the additional asymptotic term in Eq. ( pT|) to a future investigation. 

It is seen from Table [IX| that the 3/2 + state for the Gogny potential is relatively narrow. 



Taking the enhancement factors Aj = (1.0,5.3,4.6) results in a resonance position E : 



3/2^ 
3ii 



(2.88 — z0.23) MeV. For such a narrow resonance the CSM works properly for not too large 
values of the a angle (see Table [X]). Increasing it, however, the resulting energy eigenvalue 
is totally wrong what can be traced back to the increasing oscillations of the NN potential 

15 



resulting from its analytical continuation to the complex r-values. 

In the last step we checked how the results depend on the NN interaction used. We 
choose a modern, realistic NN interaction as given by the Reid 93 potential. The resulting 
positions of the 3n resonances, as far as we could trace them out, are shown for two examples 
in Table KIl In spite of the fact that the particular values of the Ay's change of course, the 
resulting picture is basically the same as for the Gogny potential. We are still too far away 
from all Aj = 1 in order to definitely say that all the widths will be quite large. But this 
will likely be the case. Again to achieve converged results one has to take at least 4 terms 



in the boundary condition (21) into account. 



IV. SUMMARY 

We investigated the configuration space Faddeev equations with the aim to determine the 
location of 3-neutron resonances. Realistic nn forces have been used. The asymptotic form 
in the 3n breakup channel is known only in the sense of an asymptotic expansion. We found 
that several terms thereof have to be included in order to get stable and reliable results, which 
are independent of the matching radius. The search for the resonance positions was carried 
through by artificially enhancing 2n force components such that a 3n bound state exists. 
Then by reducing those enhancement factors we followed the paths of the energy eigenvalue 
into the second sheet adjacent to Re(E) > 0. Unfortunately it turned out that numerical 
instabilities set in before we reached the actual strengths factors 1 for all the individually 
enhanced force components. Therefore we could not reach the very final positions of the 
studied 3n resonances in the states J n = l/2 - ,3/2~, and 3/2 + . The detailed reason for 
those instabilities remains to be found out. Despite of that our results already indicate that 
the imaginary parts of the complex resonance energies will be quite large, so that one should 
not expect to see resonances in experiments. 

We also applied the CSM, which works beautifully in atomic phyics. The complex co- 
ordinates turn the exponentially increasing resonance wave function into an exponentially 
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decreasing one, like for a bound state, and the problem with the approximate boundary 
conditions is avoided. However, for realistic NN interactions the complex scaling trans- 
formation r — > r e ia increases the interaction range by a factor — !— and leads for scaling 

o j cosa 

angles a greater than about 40° to a strong oscillations of the potential. This restricts the 
application of the complex scaling method in more than 2-nucleon systems to rather narrow 
resonances. In the 2-nucleon system one can avoid the above restriction by performing the 
complex rotation outside the region of non-zero interaction. Due to the neccessary particle 
permutation in the 3-body system we were not able to use that idea in the 3-body context. 

This project deserves further investigations to locate all lowest resonances in the three 
nucleon systems based on modern NN and possibly three-nucleon forces. 
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APPENDIX A: PERMUTATION OPERATOR 



To evaluate the term on the right hand side of Eq.(|Tq) one needs 
< xyf3\P\ip >= ^2 J dx x' 2 dy'y' 2 < xyf3\P\xy'~f > 



7 {x',y') 



x'y' 



(Al) 



where the permutation operator is P = P12P23 + P13P23 and the set of discrete quantum 
numbers abbreviated by (3 or 7 is {(ls)j(\l/2)l(jl) JM\ (tl/2)TM-r}. Here l,s, and j are 
the orbital angular momentum, total spin and total angular momentum in the two-body 
subsystem 2-3, A, 1/2, and I are the corresponding quantities for particle 1, and JM denote 
the conserved 3N total angular momentum and its magnetic quantum number. In isospin 
space, t refers to particles 2 and 3 and couples with the isospin 1/2 of particle 1 to the total 
isospin T and its magnetic quantum number Mt- Following |27j one obtains for the matrix 
element of the permutation operator 



< xya\P\x'y' a >- 



du 



5(x'-x')S(y'-y>) ~ ca/l 
.i\l'+2 (y') x '+ 2 



X' 



with 



x' — J -x 2 + xyu + y 2 



~> 9 2 3 j_ 1 2 
V = V 16 ~ ~i XyU 4 



and 



(A2) 



(A3) 



r l+ s l („A r 2+S2 ^^iSifc 
act' 



k r 1 +r 2 =l' si+S2=y 

The purely geometrical quantity g 1 a 1 a> lk * s given by 



%& k = -v 


'ii'ss'jj'XX'iPii'i 


-1) 


fc4 




1^1+^2 






' I s 


J 




' V s' 


r 


E 


-1) L LS< 


A I 


I 




A' \ 


V > 


LS 


















L S 


J. 




L S 


J 



.3. 
4' 



1/2 1/2 n r 1/2 1/2 s 

1/2 T f / \ 1/2 S s' 



(A4) 
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\ (2r 1 )\(2r 2 y.i 



E^l 



A'! 



E*3 



Vi si r 3 \ / / r 3 fc' 

/ \0 0, 

>2 s 2 s 3 \ /A s 3 



I r 3 
s 3 A L 



\ (2 Sl )\(2s 2 )\ ~ Vo o 0/V0 o o 
n r 2 I 

< si s 2 A' > (A5) 

The superscript cart stands for cartesian and I = 21 + 1. To arrive at the kernel on the 
right-hand side of Eq. (|TSf) one has to evaluate 

"i G^(x,y,«)0 7 (£',y') 



/i < 



X'hy'^i X',y' 



We introduce polar coordinates 



(A6) 



x = pcosO' 
y' = —^-psin9' 



and changing the w to the 0' integration 



4 / >6I+ 



• • • = -4= / d0' 



i x'y' 



v/3 



one arrives at 



4 r 

< xy/3\P\i(f >= E ^ y e 



{cose') l -y{^-sine') x 



-0 7 (p,^) 



with 



G 



and 



k n+r2=/ 7 si+S2=A 7 



cos 2 6>' — \cos 2 9 — ^sin 2 9 
u(9,9 ) = — 



From Eq.(|M3|) follows that 

K Pl {9,9') 



^cos9sin9 



4 G p £ ar (9,9') 
VS (cos9') L i^sin9')^ 



(A7) 



(A8) 



(A9) 



(A10) 



(All) 



(A12) 
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APPENDIX B: THE SOLUTION OF THE ENERGY EIGENVALUE PROBLEM 
INCLUDING N-TERMS IN THE ASYMPTOTIC BOUNDARY CONDITION 



Due to numerical reasons it is advantegous to extract the exponential part from the 
reduced Faddeev amplitude 4> a (p,9) and to work instead with 

F a (p,9) = ( j )ct (p,9)/(pe ik n (Bl) 

Following ref. |28| one introduces suitable grid points in p {p iy i — 0, • • • , /; p = 0, pj = 
Pmax} and in 9 {9k, k — 0, • • • , K; 9 = 0, 9k = f }• Then one can write F a (p, 9) in terms 
of cubic splines s m (p) and s n {9) as 

A/+1 N+l 

F a (p,9)= a amn s m (p)s n {9) (B2) 

m=0 n=0 

with M = 21 and iV = IK. 

Due to the properties of the cubic splines as given in |2?| one has 



F a (pi, 9k) — a a ,2i,2k 
dF Q 



dp 



A — a a,2i+l,2k 



09 

The boundary conditions at the origin (p0|) imply that 



dF a 

' Pifik = a a,2i,2k+l (B3) 



0-a,m,n=N = 0. (B4) 



Writing the asymptotic expansion (|21| ) for F a (p, 9) under the assumption that no two-body 



bound state exists in any channel a results in 



when N t terms are included in the asymptotic form. 
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Imposing the boundary condition (|B5|) at p = p max implies the dependence of a a ,m=M+i,n 
on a a ^ m= M, n and Bf(9). In order to remove the unknown Bf{9) functions the boundary 
condition has to be imposed on the last Nt points of the p-grid. This leads to the dependence 
of 

o-a,m=M+i,n on Cb a ,m=M,n-, ' • ~i o>a,m=M-2(N t -i),n- In the following the form of this dependence 
is derived. 

The asymtotic form for the partial derivative is according to (|B5| ) 



(>K:{ ' K ° ] - -A x {p)F a {p,0) -f:Mp)B?(6) (B6) 



9 P 1=2 

with 



MP) = f (B7) 
2p 



and for / > 2 



Taking the partial derivative of ( |B2|) with respect to p and the form (p6|) yields for the last 
N t points of the p grid {p = pi-i+i, i = l,---,N t } 

a a , M +l-2(i-l),nS n (0) = -^ a «,M-2(i-l),nS„(6 l )y4i(p7_i + i) - ^ Ai{p I ^ i+l )Bf{9). (B9) 
n n 1=2 

Eq. (|B"9"|) together with its ^-derivative taken at the points of the #-grid lead to 

a a ,M+l-2(i-l),n = — <VM-2(i-l),n^4l(p/-i+l) ~~ A(pr-m)-ffi*(flfc) (BIO) 

1=2 



j Da 

where for n = 2k Bf = Bf(9 k ) and for n = 2k + 1 B? = 



The set ( |B10| ) are N t equations for N t — 1 unknowns Bf{6k). They allow to write a at M+i,n 
in terms of the a a ,M-i+i,n, i — 1, • • • , 2(iV t — 1) + 1 as 

2(jV t -l)+l 

Q (Bll) 

1=1 

where the Cj's are functions of all the Ai(pj^ i+ iYs resulting from ( |B10| ). 

In this way, incorporating the boundary conditions flB"4] ) and ( P5| ) one can write F a (p, 9) 
on the chosen grids as 
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M N N 2(jV t -l)+l 

Fa(p,6)= EE (pKW + E( e 

Q>a,M— i+l,n Ci\SM+i{p)s n {9) (B12) 

m=l n=l n=l i=l 

where, due to the last equation in ( |B4| ) we renumbered n = N+ l by n = N. 

Inserting the expansion ( |B12| ) with the help of ( |B~1| ) into Eq.(|T8|) and choosing collocation 



points p p , p = 1, • • • , 21 and 9 q , q = 1, • • • , 2K |28| yields 

N c M N N c M N 

^ ] ^ ' ^ ] ^-apq,f3mn^(3mn ^ ] ^ ] ^ ] B a pq^ mn dp rnn (B13) 
/3=1 m=l n=l /3=1 m=l n=l 



with 



// 3 / 1 // 

A apq ,i3mn = 5 a p[s'^{p p )s n {6 q ) + ( h 2ik)s m (p p )s n (9 q ) + — s m (p p )s^(6> 5 ) 

Pp Pp 

f 1(1 + 1) A(A + 1) , 1 3iA; xl , , xl ,„ 1j)N 

p;cos 2 8 q p*sin*9 q p 2 p p p 



and 



Bapq^mn ~ V a p(p p COsQ q ) S m (p p ) S n (9 q ) 

+ T, v ^(pP cosd q) s m(p P ) / K lP {9 q ,d')s n {d')de' 

7=1 "' 6, 

+^m,M-i+l (l<i<2(JV*-l)+l) { — ^a/jCi [ S A/+1 {Pp) S n{9 q ) 
3 / 1 // 

+(— + 2i/c)s' M+1 (p p )s n (0 g ) + —s M +i(p P )s n (6 q ) 
Pp Pp 

f 1(1 + 1) A(A + 1) . 1 3ifc xl . . 

2Q~ + 2 ■ 20 ~(~ + )} S M+l{pp)Sn{0q)\ 

pp\cos l 9 q (f v s%n l Q q p 2 p pp 

V a/3 (ppCOsOg) S M +1 (Pp) s n {9 q ) 

+ Y. V ^PvCos6 q )s M+l { P p) / K lP {9 qi 9')s n {9')d9'} (B15) 

7=1 Je 

The number of channels is denoted by N c . This eigenvalue problem (|B13j ) can be treated in 
a standard way by changing it into the form 

A~ x Ba = Xa (B16) 

and searching for that energy at which A = 1. 

To solve Eq. (|B16| ) we used a Lanczos type iteration method |RJ. For the application 



of B and A 1 on a given vector a we applied the tensor product method . Hereby A is 
decomposed as 
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A= A® B+C®D 



with 

A P m = s m {p p )(— + ) + (— + 2ik)s' m {p p ) + s' m (p p ) 

Pp Pp Pp 

= s apSn (e q ) 

Cpm = ^''5m(Pp) 
Pp 

Then we use the identity 

A® B + C <g> D = (AU- 1 <g> <g> S + I <g> /)([/ (8) y) 

with 

iT 1 /) = V~ l EV 

and II, S diagonal. 
It follows 

(A® B + C <8) D)- 1 = (U- 1 <g> V" 1 )^ <g> E + / <g> J)" 1 ^" 1 <g> l/S- 
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TABLES 



NN state 


Ay 


Tmax 


E™ s (MeV) 








exact (m=l) 


m=0 


m=l 


m=2 




5.50 


10.0 
15.0 
20.0 


1.42-i0.75 
1.42-i0.75 
1.42-i0.75 


1.75-il.04 
1.58-i0.63 
1.36-i0.58 








5.00 


10.0 
15.0 
20.0 


15.08-i0.80 
15.08-i0.80 
15.08-i0.80 


15.11-i0.90 
15.07-i0.85 
15.06-i0.83 


15.08-i0.80 
15.08-i0.80 
15.08-i0.80 




3 P 2 - 3 F 2 


6.20 


10.0 
15.0 
20.0 


2.79-il.l7 
2.79-il.l7 
2.79-il.l7 


3.07- il.05 
2.62-i0.99 

3.08- il.28 


2.79-il.l7 
2.79-H.17 
2.79-il.l7 





TABLE I. Resonance positions in different NN states obtained with the correct (column 4) 
and approximate (columns 5-7) boundary conditions applied at different matching radii r max - The 
Reid93 NN potential multiplied with the factors Ay has been used. 
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r-matching 


Re(a) 


Im(a) 


Z7127V 

res 


(fin) 


(°) 


(°) 


(MeV) 


10.0 


0.0 


0.0 


1.75 - i 1.04 


10.0 


20.0 


0.0 


1 40 - i 88 


10.0 


0.0 


-40.0 


1.36 - i 0.58 


10.0 


0.0 


-80.0 


1.24 - i 0.58 


10.0 


20.0 


-40.0 


1.42 - i 0.73 


10.0 


20.0 


-80.0 


1.42 - i 0.75 


2.0 


20.0 


-80.0 


1.22 - i 0.93 


2.0 


20.0 


-160.0 


1.43 - i 0.75 


20.0 


20.0 


0.0 


1.42 - i 0.73 


20.0 


30.0 


0.0 


1.43 - i 0.74 


20.0 


40.0 


0.0 


1.43 - i 0.75 


20.0 


44.0 


0.0 


1.48 - i 0.75 


20.0 


46.0 


0.0 


1.71 - i 0.72 


20.0 


48.0 


0.0 


2.99 - i 0.36 



TABLE II. Complex energy eigenvalues for the 3 Po 2N state obtained with the Reid 93 NN 
interaction enhanced by Ay = 5.5 and an approximate boundary condition (m = term only in 
Eq.(|8|) ). The correct energy eigenvalue is = (1.42 — i 0.75) MeV. 
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r-matching 


A 


Tn 
1 u 


a 


E 2N 

^res 


(fin) 


(fm~l) 


( fm) 


(°) 


(MeV) 


10.0 


2.0 


5.0 


10.0 


15.08 - i 0.84 


10.0 


2.0 


5.0 


20.0 


15 07 - i 81 


10.0 


2.0 


5.0 


30.0 


15.07 - i 0.80 


10.0 


2.0 


5.0 


40.0 


15 08 - i 80 


10.0 


2.0 


5.0 


50.0 


15 08 - i 80 


10.0 


2.0 


5.0 


60.0 


15.08 - i 0.80 


10.0 


2.0 


5.0 


70.0 


15 08 - i 80 

A '.V • VA v. / l v / • v. / v / 


10.0 


2.0 


5.0 


80.0 


15.08 - i 0.80 


15.0 


2.0 


5.0 


60.0 


15.08 - i 0.80 


20.0 


2.0 


5.0 


60.0 


15.08 - i 0.80 


15.0 


2.0 


3.0 


60.0 


15.08 - i 0.80 


15.0 


3.0 


3.0 


60.0 


15.08 - i 0.80 



TABLE III. Smooth-exterior complex scaling applied to the Do, resonance state of the Reid93 
NN potential with Ay = 5.0. Only the first term in the boundary condition @ was taken into 
account. The correct energy eigenvalue is E^es = (15.08 — i 0.80) MeV. 
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NN states 


N r 
* C 


3 H binding energy (MeV) 






p-space 


r-space 


% + 3 Si - 3 D 1 


5 


-7.630 


-7.631 


i<i 


10 


-7.403 


-7.406 


3 < 2 


18 


-7.747 


-7.750 



TABLE IV. Triton binding energies obtained in momentum- and coordinate-space calculations 
with the Reid 93 NN potential using different numbers of channels iV c . 
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/9-matching 


Re(a) 


Im(a) 


rp3B 

res 


ffm) 


(°) 


(°) 


(MeV) 


10.0 


0.0 


0.0 


8.50 - i 1.81 


10.0 


20.0 


0.0 


8.48 - i 1.87 


10.0 


40.0 


0.0 


8.46 - i 1.86 


5.0 


20.0 


0.0 


8.36 - i 1.78 


5.0 


40.0 


0.0 


8.71 - i 2.03 


5.0 


40.0 


-20.0 


8.45 - i 1.93 


5.0 


40.0 


-40.0 


8.46 - i 1.86 


5.0 


40.0 


-60.0 


8.47 - i 1.86 



TABLE V. CSM applied to a 3-boson system interacting through the S-wave interaction of 
Eq.dH) in the J 71 " = 0+ state. The correct energy eigenvalue is Eff s = (8.47 - i 1.86) MeV. The 



boundary condition (21) with A and B terms only has been applied. 
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Q>mnv (fill) 


E^f s (MeV) for particular number of terms in Eq.(21 


) 




1 


2 


3 


4 


10.0 


8.40-il.40 


8.50-11.81 


8.48-H.85 


8.49-il.85 


15.0 


8.81-il.80 


8.46-il.87 


8.46-H.85 


8.46-il.85 


20.0 


8.00-i2.18 


8.46-il.85 


8.46-H.85 


8.46-il.85 


30.0 


8.64-il.69 


8.47-il.86 


8.47-H.86 


8.47-il.86 



TABLE VI. Resonance positions of a 3-boson system in the J w = + state interacting through 
the S-wave interaction of Eq.(23). The results are without complex scaling but taking an increasing 
number of terms in the boundary condition (^) into account. The correct energy eigenvalue is 
Efes = (8-47 — i 1.86) MeV. Number of terms 1 means that only the A{6) term has been used. 
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E 2r , (MeV) 


E-*„ (MeV) 


1 n 




3 P 2 ~ ^2 






1.8360 


9.2804 


5.2012 


-4.0 


-31.68 


1.7241 


8.9456 


5.1422 


-3.0 


-28.06 


1.5960 


8.5851 


5.0797 


-2.0 


-24.15 


1.4370 


8.1846 


5.0123 


-1.0 


-19.73 


1.3297 


7.9584 


4.9758 


-0.5 


-17.09 


1.1930 


7.7518 


4.9442 


-0.1 


-14.21 


1.1205 


7.6981 


4.9365 


-0.01 


-12.96 


1.05 


7.65 


4.90 




-11.13 


1.0 


7.5 


4.8 




-8.24 


1 n 

J. .u 


1 n 

i-.VJ 


4 8 




-7.12 


1.0 


1.0 


4.7 




-5.13 


1.0 


1.0 


4.6 




-3.29 


1.0 


1.0 


4.5 




-1.65 


1.0 


1.0 


4.4 




-0.23 



TABLE VII. Required strength factors Ay for particular force components of the Gogny NN 
interaction together with the corresponding 2- and 3-neutron (3/2~) bound energy eigenvalues. 
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/VJ P2— J i<2 


E 3n (Me V) for particular number of terms in Eq. (^T|) 




1 


2 


3 


4 


5 


6 


4.30 


0.86 - i 0.068 


0.86 - i 0.068 


0.85 - i 0.092 


0.85 - i 0.092 


0.86 - i 0.068 


0.86 - i 0.068 


4.25 


1.26 - i 0.27 


1.34 - i 0.21 


1.33 - i 0.21 


1.34 - i 0.21 


1.33 - i 0.22 


1.33 - i 0.22 


4.20 


1.88 - i 0.38 


1.77 - i 0.41 


1.77 - i 0.41 


1.77 - i 0.41 


1.79 - i 0.38 


1.79 - i 0.38 


4.15 


2.24 - i 0.53 


2.20 - i 0.62 


2.19 - i 0.62 


2.19 - i 0.63 


2.21 - i 0.64 


2.19 - i 0.63 


4.10 


2.54 - i 0.70 


2.58 - i 0.87 


2.60 - i 0.86 


2.58 - i 0.88 


2.63 - i 0.89 


2.57 - i 0.84 


4.05 


2.72 - i 0.90 


2.94 - i 1.15 


2.96 - i 1.16 


2.97 - i 1.14 


3.04 - i 1.15 


2.92 - i 1.18 


4.00 


2.81 - i 1.08 


3.27 - i 1.49 


3.32 - i 1.46 


3.32 - i 1.46 


3.31 - i 1.47 


3.31 - i 1.47 


3.95 


2.84 - i 1.20 


3.63 - i 1.86 


3.64 - i 1.81 


3.66 - i 1.79 


3.64 - i 1.80 


3.64 - i 1.80 


3.90 


2.84 - i 1.30 


4.01 - i 2.24 


3.97 - i 2.16 


3.95 - i 2.16 


3.96 - i 2.14 


3.97 - i 2.14 


3.85 


2.84 - i 1.37 


4.39 - i 2.55 


4.25 - i 2.53 


4.24 - i 2.53 


4.22 - i 2.54 


4.21 - i 2.51 


3.80 


2.84 - i 1.41 


4.70 - i 2.82 


4.54 - i 2.91 


4.50 - i 2.92 


4.47 - i 2.94 


4.49 - i 2.93 


3.75 


2.84 - i 1.44 


4.88 - i 3.13 


4.77 - i 3.29 


4.76 - i 3.32 


4.73 - i 3.33 


4.78 - i 3.35 


3.70 


2.84 - i 1.47 


5.07 - i 3.40 


4.99 - i 3.68 


4.95 - i 3.75 


4.98 - i 3.78 


4.96 - i 3.77 


3.65 


2.84 - i 1.50 


5.17 - i 3.66 


5.12 - i 4.08 


5.14 - i 4.18 


5.16 - i 4.17 


5.14 - i 4.19 


3.60 


2.84 - i 1.51 


5.23 - i 3.91 


5.21 - i 4.51 


5.30 - i 4.67 


5.31 - i 4.68 


5.33 - i 4.69 



TABLE VIII. The 3-neutron 3/2 _ resonance state trajectories for the Gogny NN interaction 
obtained with increasing number of terms in the boundary condition (f2~l|). The enhancement factors 
A for the 1 S'o and 3 Pq nn force components are equal one. 
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J 77 


Ay 


E-Kr, fMeV) 




So 








1/2- 


1.0 


5.5 


2.75 


2.82 - i 2.82 


3/2" 


1.0 


1.0 


3.70 


4.95 - i 3.75 


3/2+ 


1.0 


5.0 


4.3 


5.74 - i 1.53 



TABLE IX. The positions of 3-neutron resonances for the Gogny potential obtained with four 
leading terms in the boundary condition (f2l|). 
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Re(a) 


Im(a) 


'Oil 


(°) 


(°) 


(MeV) 


10.0 


0.0 


2.88 - i 0.23 


20.0 


0.0 


2.88 - i 0.23 


30.0 


0.0 


3.11 - i 0.16 



TABLE X. CSM applied to the J w = 3/2 + resonance state of the Gogny potential with en- 

3/2+ 

hancement factors A« = (1.0, 5.3, 4.6). The correct energy eigenvalue is E 3n = (2.88— i 0.23) MeV. 
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Ay 


E 3n (MeV) 




1 S'o 








3/2" 


1.0 


1.0 


3.25 


5.30 - i 3.53 


3/2+ 


1.0 


4.6 


3.50 


5.93 - i 1.55 



TABLE XL The positions of 3-neutron resonances for the Reid 93 potential obtained with 



four leading terms in the boundary condition ( |2l| ) 
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FIGURES 

FIG. 1. Analytical continuation of the 3 Po Reid93 NN potential under the complex scaling 
transformation r —* r e' ia for five values of a = 0°, 20°, 40°, 60° and 80°. ((a) real part, (b) 
imaginary part). 

FIG. 2. The real and imaginary parts of the reduced Faddeev amplitudes F for a 3-boson 
resonance in the state + without complex scaling ((a) real part, (b) imaginary part). 

FIG. 3. The same as in Figj| but with ECSM taking a = (40.0° - z40.0°) ((a) real part, (b) 
imaginary part). 

FIG. 4. The 3-neutron 1/2 - , 3/2~, and 3/2 + resonance state trajectories for the Gogny NN 
interaction obtained with 4 terms in the asymptotic boundary condition. Different points corre- 
spond to resonance positions with different values of the enhancement factors. For the l/2~ state 
Ai£ = 1.0, \3p changes from 7.5 to 5.5 in steps of 0.1, and A3p 2 _3p 2 changes from 4.75 to 2.75 
in steps of 0.1. For the 3/2~ state \i So = 1.0, A3p = 1.0, and A3p 2 _ 3 p 2 changes from 4.3 to 3.7 
in steps of 0.05. For the 3/2 + state \i So = 1.0, A3p changes from 5.5 to 5.0 in steps of 0.05, and 
A3p 2 _3 F2 changes from 4.8 to 4.3 in steps of 0.05. 
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